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The breakup of shearless invariant tori with winding number uj — [0, 1, 11, 1,1,.. .] (in continued 
fraction representation) of the standard nontwist map is studied numerically using Greene's residue 
criterion. Tori of this winding number can assume the shape of meanders (folded-over invariant tori 
which are not graphs over the x-axis in (x, y) phase space), whose breakup is the first point of focus 
■ here. Secondly, multiple shearless orbits of this winding number can exist, leading to a new type 

' of breakup scenario. Results are discussed within the framework of the renormalization group for 

| area-preserving maps. Regularity of the critical tori is also investigated. 

c : 

In recent years nontwist maps, area-preserving maps that violate the twist condition locally in phase 
\ space, have been the subject of several studies in physics and mathematics. These maps appear 
naturally in a variety of physical models. An important problem is the understanding of the breakup 
of invariant tori, which correspond to transport barriers in the physical model. We conduct a detailed 
study of the breakup of two types of invariant tori that have not been analyzed before. 
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I. INTRODUCTION 

We consider the standard nontwist map (SNM) M as introduced in Ref. [J 



x n +i = x n + a(l- yl +1 ) 



where (x,y) € TxR are phase space coordinates and a, b £ R are parameters. This map is area-preserving and 



Un+i = y n -bsin(2irx n ) , (1) 

o 

, violates the twist condition, dx n j r \{x n ,y n ) / dy n ^ 0, along a curve in phase space. Although the SNM is not generic 
due to its symmetries, it captures the essential features of nontwist systems with a local, approximately quadratic 
\ extremum of the winding number profile. 
\& • Nontwist maps are low-dimensional models of many physical systems, as described in Refs. 121314. Of particular 
interest from a physics perspective is the breakup of invariant tori (which we alternatively call invariant curves), 
consisting of quasiperiodic orbits with irrational winding number (see Appendix) , that often correspond to transport 
barriers in the physical system. 

One important characteristic of nontwist maps is the existence of multiple orbit chains of the same winding number. 
For the SNM, in particular, the symmetry S(x, y) = (x + 1/2, —y) guarantees that whenever an orbit chain of a certain 
winding number exists, a second chain with the same winding number can be found. 

Changing the map parameters a and b causes bifurcations of periodic orbit chains with the same winding number. 



5h I Orbits can undergo stochastic layer reconnection ( "separatrix" reconnection), or they can collide and annihilate. The 
simplest reconnection-collision scenarios, which involve sequences of collisions between elliptic and hyperbolic orbits 
of a single pair of periodic orbit chains with the same winding number, have been known since early study of nontwist 
systems^ In the SNM, two standard scenarios can be distinguished that describe, respectively, the reconnection- 
collision sequence for a pair of either even or odd-period orbit chains. A detailed review of these scenarios, as well as 
a discussion of earlier studies of reconnection-collision phenomena in theory and experiments, can be found in Ref. Q. 

As recently reported^ even the simple, non-generic SNM can have more than two orbit chains of the same winding 
number, and thus reconnection-collision scenarios are more intricate than previously expected. To illustrate this, two 
stages of the odd-period standard scenario are shown in Fig. ^ Each of the winding number profiles exhibits two 
peaks at or somewhat below the winding number of the colliding odd-period orbit, which is marked by the dashed 
vertical line. Between the peaks lies a recess, i.e., for any winding number between the lowest value of the recess and 
the peak values, e.g., the value marked by dotted verticals, four (or more) orbits of that winding number traverse the 
x = 0.5 symmetry linei^ Similar observations can be made for other symmetry lines. Therefore, whenever a periodic 
orbit is studied with winding number inside the recess associated with a major odd-period orbit collision, complicated 
reconnection-collision scenarios are possible, which were called non-standard scenarios in Ref. 0- 



FIG. 1: Phase space (left) and winding number profile along x — 0.5 (right) of two stages of odd-period standard scenario. 
The y-ranges of all plots are identical. Parameters were chosen (a) slightly before collision and (b) slightly after collision of two 
major odd-period orbit chains (winding number marked by dashed lines). 



When two quasiperiodic orbits collide, the winding number profile shows a local extremum and the orbit at collision 
is referred to as the shearless curve. In previous studies of the SNM only the shearless curves invariant under the 
full symmetry group Q of the SNM (composed of the symmetry S as well as the involutions I\ and have been 

considered. However, as seen in Fig. ^b), in addition to the central extremum (here a minimum) in the winding 
number profile, other shearless curves (here marked by the two outer peaks) may exist. Figure D^a) also exhibits 
peaks, but since the plateau and spike are associated with an elliptic point and the invariant manifolds of hyperbolic 
points, respectively, rather than quasiperiodic orbits, they do not qualify as shearless curves. From now on, we will 
refer to the C/-invariant curve as the "central shearless curve" and to others as "outer shearless curves." Since the 
breakup of outer shearless curves has not been studied so far, this will be the focus of our investigation in Sec. IIVI 

Another consequence of the violation of the twist condition is the occurrence of meanders, quasiperiodic orbits that 
arc "folded over", i.e., not graphs y(x). Whereas Birkhoff's theorem states that such curves cannot exist in twist 
maps, they can occur in nontwist maps. In the SNM, meanders appear between (in parameter space) the reconnection 
and collision of odd-period orbits^ Figure ^a) shows an example, whereas in Fig. ^b), the meander has changed to 
a graph again. As seen in Fig. [Ja), the region in which meanders are found corresponds to a recess in the winding 
number profile. However, the converse is not true: Fig.^b) shows an example where meanders are absent, but still 
a recess in the winding number profile is observed. To our knowledge, the breakup of meanders has not been studied 
previously. 

The paper is organized as follows. In Sec.[n]we review previous results about the SNM relevant to this investigation. 
Section III Al contains a brief discussion of the parameter space of the SNM, in particular the details of where the 
scenarios studied in this paper occur. Section III Bl contains an account of how Greene's residue criterion is used for 
detecting the breakup of critical invariant tori. In Sec. lIIII we discuss the results for the breakup of the central shearless 
meander of winding number tu = [0, 1, 11, 1, 1, . . .] (in continued fraction representation), while in Sec. UVI we consider 
the breakup of the outer shearless tori of the same winding number. Questions of regularity of these critical tori and 
a comparison with previous results are addressed in Sec.0 Finally, Sec. IVII contains our conclusion and a discussion 
of open questions. Basic definitions are given in the Appendix. 

II. BACKGROUND 
A. Parameter space overview 

In order to identify the parameter regions where the bifurcations described in Sec.[I]occur, various parameter space 
curves can be computed (usually numerically). Collisions of periodic orbits are described by the bifurcation curves^ 
introduced in Ref. |^ and generalized in Ref. |2j. Reconnections do not occur precisely on parameter space curves (see, 
e.g., Ref.0), but within a finite range of parameters; however, the range is usually small enough that the method of 
Ref. |S| (implemented in Ref. 0) yields curves that represent a good approximation of the reconnection thresholds for 
odd-period orbits. For even-period orbits, reconnection coincides with the collision of hyperbolic orbits. 

By numerically computing the branching points at which bifurcation curves for various higher periodic orbits split 
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FIG. 2: Parameter space overview of SNM, with higher orbit branching points, bifurcation (collision) curves, and approximate 
reconnection thresholds for odd-period orbits 2/3 and 1/1. Bifurcation curves for even-period orbits 1/2, 5/6, and 11/12 and 
the breakup boundary for the central shearless tori are indicated. New points are marked by two circles. (0 denotes outer 
shearless curve, all others are central ones.) 



up into several branches below a major odd-period orbit collision, one obtains a good estimate of the parameter region 
for which multiple orbit chains exist. For an extensive discussion of these curves and their computation see Ref. |2- 

An overview of parameter space with thresholds for several examples of low-period orbits is shown in Fig. [21 
Higher orbit branching points, bifurcation curves, and approximate reconnection thresholds are shown for the odd- 
period orbits 2/3 and 1/1. Meanders occur between reconnection and collision; a recess in winding number profile 
is encountered in the region limited by branching points and bifurcation curve. For even-period orbits with winding 
number 1/2, 5/6, and 11/12, only bifurcation curves are shown, since even-period orbits do not induce branching and 
reconnections coincide with collisions. Of the corresponding winding numbers, no orbits exist above the highest (in 
6) bifurcation curve and two orbit chains arc found below the lowest one; in between various numbers of orbits exist 
on each symmetry line. In addition, the figure also contains the ragged breakup boundary, introduced in Ref. 0, 
above which the central shearless orbits have become chaotic^ We further indicate po ints (by triangles) for which the 
breakup of the central shearless curve has been studied in detail in the past (see Refs. l3l4llllTll2lll3l) as well as the two 
points investigated in this paper: The central meander from Sec. Illll is shown as a solid circle (•), the outer shearless 
curves from Sec. IIVI as an empty circle (©). Note that all breakup points for central shearless tori are located on the 
breakup boundary, whereas the outer shearless curves break up at smaller parameter values. 

The winding number investigated here was chosen such that the central shearless curve is a meander, i.e., the point 
• is located between the 1/1 reconnection and collision thresholds, and that multiple orbit chains and hence multiple 
shearless curves can occur, i.e., the point is located to the right of the 1/1 branching threshold. This can be ensured 
by picking a winding number close to a periodic orbit, here 11/12, whose bifurcation curve both branches due to a 
nearby major odd-period orbit, and has one branch crossing a collision threshold of this major odd-period orbit before 
crossing the breakup boundary. 

B. Greene's residue criterion 

Whereas the breakup boundary in Fig.[2|provides a rough estimate of the parameter values at which central shearless 
tori break, a significantly more precise tool for studying the breakup of a particular torus with given winding number 
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is provided by Greene's residue criterion, originally introduced in the context of twist maps^ This method relics on 
the numerical observation that the breakup of an invariant torus with irrational winding number lo is determined by 
the stability of nearby periodic orbits. Some aspects of the validity of this criterion have been proved for nontwist 
mapsii^ 

To study the breakup, one considers a sequence of periodic orbits with winding numbers m^/n^ converging to u>, 
limj—too mi/rii = lo. The sequence converging the fastest, and hence the most commonly used one, consists of the 
convergents of the continued fraction expansion of lo, i.e., mi/rii = [ao, ai, . . . , a*], where 

w = [oo,ai,a2, . . •] = ao H • (2) 



The stability of the corresponding orbits can be expressed through their residues, R4 = [2 — Tr(Z?M ni )]/4, where Tr 
is the trace and DM ni is the linearization of the n, times iterated map about the periodic orbit: An orbit is elliptic 
for < R4 < 1, parabolic for Ri = and R t — 1, and hyperbolic otherwise. The convergence or divergence of the 
residue sequence associated with the chosen periodic orbit sequence then determines whether the torus exists or not, 
respectively: 

• lim^oo = if the torus of winding number lo exists. 

• lim^oo =00 if the torus of winding number lo is destroyed. 

At the breakup itself, various scenarios can be encountered, depending on the class of maps and invariant torus under 
consideration. 

For twist systems, this criterion has been used to study the breakup of noble invariant tori in the standard (twist ) 
map, i.e., orbits with winding numbers that have a continued fraction expansion tail of l's (see, e.g., Refs. llJldll7h . 
It was found that at the point of breakup the sequence of residues converges to either R^ m 0.25 or a three-cycle 
containing 0.25 ... as one of its elements. 

In the standard nontwist map, the residue criterion was first used in Ref. to study the breakup of the central 
shearless torus of inverse golden mean 1/7= (VE — l)/2 = [0, 1, 1, 1, . . .] winding number. There it was discovered 
that the residue sequence converges to a six-cycled Similar studies were conducted for noble central shearless tori of 
winding numbers lo = I/7 2 (Refs. I MTT|) and lo = [0, 2, 2, 1, 1,1,.. .] (Rcf. IT^l . and the same six-cycle was found. The 
parameter values at which these shearless tori break, i.e., at which six-cycles of residues are encountered, are marked 
by triangles in Fig. [21 In this paper, we study the noble winding number w = [0,1,11,1,1,1,...], where the large 
number 11 in the second convergent had to be chosen to ensure that the breakup occurs in a region in parameter 
space where both meanders and multiple shearless tori are possible, as described in Sec. Ill AI 

In addition to nontrivial residue convergence behavior, invariant tori at breakup exhibit scale invariance under 
specific phase space re-scalings. All these results suggest that certain characteristics of the breakup of noble invariant 
tori are universal, i.e., the same within a large class of area-preserving maps. To inter pret t he results, a renormalization 
group framework based on the residue criterion has been developed (see, e.g., Refs.aiaUJUJ)- 

III. BREAKUP OF THE lu = [0, 1, 11, 1,1, .. .] CENTRAL SHEARLESS MEANDER 
A. Search for critical parameter values 

In order to study a shearless irrational orbit, one needs to locate parameter values on its bifurcation curve. 
This can be achieved numerically by approximating them by parameter values on the bifurcation curves of nearby 
periodic orbits, usually of orbits with winding numbers that are the continued fraction convergents of lo. For 
lo = [0, 1, 11, 1,1,...] w 0.920748351059159504, the convergents up to the highest numerically accessible one in our 
studies are shown in Table [I] 

For given parameters (a, 6), any of these periodic orbits (if they exist) can be found along symmetry lines via a 
one-dimensional root search, as explained, e.g., in Ref. 0. Performing this search for a range of parameters, usually 
varying b while keeping a constant, results in the relation y(b), i.e., the location(s) of the periodic orbit along a given 
symmetry line, as shown in Fig.[3]for the orbits 11/12, 12/13, 23/25, and 35/38 along the s\ symmetry line. 

In this plot, orbit collisions are found where two y(b) branches meet, i.e., at extrema of b(y). Especially for higher 
period orbits, multiple collisions can be observed, but in this section we focus only on the ones approximating the 
central shearless curve, deferring the primary outer ones (i.e., the only additional collisions found for the lowest, 11/12, 
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[n] 





1/1 


13 


2707/2940 


26 


1410348/1531741 


1 


11/12 


14 


4380/4757 


27 


2281991/2478409 


2 


12/13 


15 


7087/7697 


28 


3692339/4010150 


3 


23/25 


16 


11467/12454 


29 


5974330/6488559 


4 


35/38 


17 


18554/20151 


30 


9666669/10498709 


5 


58/63 


18 


30021/32 605 


31 


15640999/16987268 


6 


93/101 


19 


48575/52756 


32 


25307668/27485977 


7 


151/164 


20 


78596/85361 


33 


40948667/44473245 


8 


244/2 65 


21 


127171/138117 


34 


66256335/71959222 


9 


395/429 


22 


205767/223478 


35 


107205002/116432467 


10 


639/694 


23 


332938/361595 


36 


173461337/188391689 


11 


1034/1123 


24 


538705/585073 


37 


280666339/304824156 


12 


1673/1817 


25 


871643/946668 


38 


454127676/493215845 



TABLE I: Continued fraction convergents for [0, 1, 11, 1, 1, . . .], where [n] = [0, ao, . . . , 0,1+2] (following the notation of Ref. 0). 




b 



FIG. 3: j/-values in phase space of the 11/12, 12/13, 23/25, and 35/38 central orbits on the si symmetry line as function of b, 
at a = 1.0645342893. 



orbit) to Sec. IIVI In contrast to previous publications, in which central collisions appear as maxima in b(y), here, in 
the meander regime, they are associated with minima. 

The parameter values 6r„i (a) of these central collisions, found by extremum searches and marked by solid circles (•) 
in Fig- [31 converge to b Qa (a) (located on the bifurcation curve of ui). Now Greene's residue criterion, as described in 
Sec. Ill Bl can be used to determine whether at (a, b 00 (a)), the shearless curve still exists or not: At parameter values a 
and the best known approximation to 600(a), the residues of all periodic orbits of convergents that have not collided, 
here the orbits [n] with even n, are computed. Their limiting behavior for n — > 00 reveals the status of the torus. 
By repeating the procedure for various values of a, with alternating residue convergence to and 00, the parameter 
values of the shearless torus breakup, (a c , 6oo(a c )), can be determined to high precision. 

Due to numerical limitations, the highest orbit collision used here to approximate 600(a) is 6(33] (a). However, a 
better approximation can be obtained by observing (in hindsight) that close to the critical breakup value, the 6[„j (a) 
obey a scaling law 

6[„] = boo + B(n)5- n , (3) 

where B(n) is empirically found to be periodic in n with period 12 as n — > 00. As 630 is unknown, the scaling is most 
readily observed by plotting 

m (V+i] ~ Vl) = - « ln(5 i > ( 4 ) 

where B(n) = \n(B(n + 1)/Si — B(n)) is also periodic in n with period 12. This is shown in Fig. @] where for clarity 
only the offsets of In (6[„ +1 ] — 6r n i) about the average slope are shown. The 6[„] values used here were obtained from 
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FIG. 4: 12-cycle of 6[„] differences in approximating the critical shearless [0, 1, 11, 1,1,.. .] meander at a c — 1.0645342893 on 
si. 



orbits colliding on the si symmetry line, although the same behavior is observed on the S2 symmetry line. For S3 
and Si, a similar plot is found, however, the 12-cycles are shifted by n ± 6. The slope was calculated from the last 24 
difference values by averaging the last 12 slopes (hi(6[„ +13 ] — 6[ n +i2]) — l n (^[n+i] — ^[n])) /12, with n = 8, . . . , 20. The 
result is log<5i w 0.98496 ± 0.00036, or 

S 1 = 2.678 ±0.001 . (5) 

The periodicity of B(n) enables us to obtain a better approximation, b* (i.e., closer to than 6[ 33 j) from lower 
6[„]-valucs, using the extrapolation 

b * = b m + 77 , b[32 r?r — 1—) x « ( 6 ) 

(,0[21] ~ 0[20]J ~ (°[33] - °[32]J 

for the best shearless 6 00 (a)-approximation at which to apply Greene's residue criterion. 



B. Residue six-cycle at breakup 



Searching along (a,6*(a)) for the transition between residue convergence to and 00, we obtain as the critical 
parameters for the shearless meander breakup: 

[a c , b c ) = (1.0645342893, 0.209408148327230359) . (7) 

At these parameters, only orbits with n even in Table [I] exist, two for each n, denoted as "up" and "down" for larger 
and smaller y-values on a symmetry line, respectively. Plotting the residues of these orbits, one observes the 6-cycles 
in Fig. here shown for the up and down orbits on si. 

The same cycles are found for the other symmetry lines, with up and down orbits interchanged and shifted by 
[n] — > [n + 6] , as summarized in Table [H] Since these cycles are the same as the ones observed in Ref. (up to an 
interchange of up and down orbits and a shift of [n] — * [n + 7] ) , the labels Ci and Di here were assigned to reflect this 
correspondence. 

The convergence of the twelve individual values of the six-cycles, of which due to symmetry properties of the map 
only six are truly distinct as listed in Table ITT1 are shown in Fig. Table ITTT1 gives the corresponding numerical values. 
In addition to the residues at the critical parameters for the shearless meander breakup, at a c = 1.06453428930 and 
b c = 0.209408148327230359 (bold), residues slightly below, at a_ = 1.06453428925 and b- = 0.209408148282494088 
(dashed), and slightly above, at a+ = 1.06453428935 and b+ = 0.209408148371966630 (dotted) are displayed. All 
these parameters are extrapolated values from the bifurcation curves of the [20], [21], [32], and [33] orbits. For 
comparison, the thin solid line shows residue behavior without extrapolation, for a c = 1.06453428930 and 633 (a c ) = 
0.209408148327230605 at the [33] orbit collision. 

In summary, the six independent residues are 

C x = -0.6090 ±0.0046, C 2 = -1.2901 ± 0.0007, 
C 4 = 1.5945 ± 0.0022, C 5 = 2.3434 ± 0.0018, 
C 6 = 2.5919 ± 0.0023, D 6 = -2.6365 ± 0.0007, 
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[n] 

FIG. 5: Residue six-cycles for orbits si at meander breakup. Cycles for orbits on s 2 are the same, with up and down 
interchanged. Cycles on S3 are shifted by three, with up and down interchanged. Cycles on S4 are shifted by three. 



[n] 


Rul — Rdl 


Ru2 — Rdl 


Rui — RdA 


Rui — Rdi 


[8], [20], [32] 


a 


Di 




Di 


[10], [22], [34] 


c 2 


D 2 


c 5 


D 5 


[12], [24], [36] 


c 3 


D 3 


C 6 


D 6 


[2], [14], [26] 


d 


Di 


Ci 


Di 


[4], [16], [28] 


c 5 


D B 


c 2 


D 2 


[6], [18], [30] 


c 6 


D 6 


d 


D 3 



TABLE II: Period-six convergence pattern of residues near criticality on different symmetry lines (following the notation 
of Ref.0). Symmetry properties of the SNM further imply C e = C 3 , Di = Ca, D 2 = C 2 , D 4 = Ci, D 5 = C 2 , D 6 = D 3 . 

where these numerical values were calculated each as the average of the last four corresponding values in Table 1 1 1 1 1 
and the error given here is the standard error (standard deviation of the mean). 

C. Phase space scaling invariance and renormalization results 

The phase space at the critical parameter values for the meander breakup is shown in Fig. [7] As in previous studies, 
the shearless meander at breakup is scale invariant under specific re-scalings of phase space 

This is readily seen by zooming in at a certain point using different levels of magnification: For example, we zoom 
in on the intersection (x s ,y s ) of the shearless meander with the S3 symmetry line, and transform to symmetry line 
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FIG. 6: Convergence of residue values at meander breakup. Shown are the two corresponding pairs of residues at the critical 
parameter values (a c ,b c ) (bold line), at (o_,6_) (dashed line), and (a+,6+) (dotted line), all using extrapolated value b* (a) 
from the bifurcation curve ft^jvalues of the [20], [21], [32], and [33] orbits. The thin line shows residue behavior without 
extrapolation for (a c , 633(00)). 
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-2.1346 


-2.1346 


c 4 


1.1143 


1.1143 


1.1143 
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-0.1976 


-0.1976 
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1.6009 
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1.6009 


[20] 




-0.6151 


-0.6151 


-0.6151 




1.5925 


1.5925 


1.5925 


[26] 




-0.6125 
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-1.2897 


-1.2903 
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2.3389 


2.3402 


2.3414 


[34] 




-1.2886 


-1.2911 


-1.2940 




2.3303 


2.3455 


2.3612 
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-2.4643 


-2.4643 


-2.4643 


c 6 


2.3155 


2.3155 


2.3155 


[12] 




-2.5402 


-2.5402 


-2.5402 




2.7032 


2.7032 


2.7032 


[18] 




-2.6371 


-2.6371 


-2.6372 




2.5938 


2.5939 


2.5939 


[24] 




-2.6368 


-2.6370 


-2.6372 




2.5935 


2.5937 


2.5939 


[30] 




-2.6346 


-2.6374 


-2.6402 




2.5921 


2.5950 


2.5980 


[36] 




-2.6434 


-2.6346 


-2.6284 




2.6148 


2.5850 


2.5589 



TABLE III: Numerical values of the six independent residues. 
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X 

FIG. 7: Phase space at the [0, 1, 11, 1, 1, 1 . . .] meander breakup, a c = 1.0645342893 and b c = 0.209408148327230359. Also 
shown are the symmetry lines si, S2, S3, and S4. 



coordinates, in which the S3 symmetry line becomes a straight line, 



x' = x--(l- y z ) 

y' = y- y s ■ 



(8) 



Here, y s = 0.59790858154 was obtained by applying the same scaling as in Eq. ijBJ) for 6r n i values to the yr n i locations 
of the [22], [24], [34], and [36] periodic orbits on S3, once for the respective up orbits and once for the down orbits, 
and averaging over both results. 

Figure [S] shows two levels of magnification of the meander in these coordinates, each along with the up and down 
periodic orbits of one of its convergents. The plotted region was chosen to allow a direct comparison with Fig. 7 of 
Ref. 0- Although the two plots deviate slightly from each other towards the edges of the plotted regions (because in 
contrast to Ref. the x' and y' ranges are larger here, i.e., scales at which the meander is still influenced mostly by 
lower convergents), they correspond exactly around the origin. 

For a quantitative analysis, we compute the scaling factors a and (3 such that the meander in the vicinity of its 
intersection with S3 is invariant under (x' , y') — * (a 12 x', f3 12 y') (following the notation of Refs.00). Again, these are 
found from the limiting behavior of convergent periodic orbits: Denoting by (x' n ±,y' n ±) the symmetry line coordinates 
of the point on the up (+) or down (— ) orbit of the [n]th convergent that is located closest to (0, 0), we compute 



b n,± 



b n+12,± 



al2 
> Pn,± 



Vn,± 



y-n+l2,± 



(9) 



Averaging the six values ajg-t, «20±> and a\\±, we find a 12 — 321.65±0.070, i.e., a = 1.61789±0. 00003 (with the error 
being the standard deviation of the mean). Similarly, for /3, we obtain /3 12 = 431.29±0.19, i.e., /? = 1.65792±0.00006. 
Thes e are the scaling factors used in Fig. |H1 Within numerical accuracy, they coincide with the values found in 
Refs.l3ll8l 

To interpret the scaling invariance of the shearless meander itself and its convergents under [n] — > [n + 12] one 
can introduce a renormalization picture, with a renorma lization group operator IZ^ acting on the space of maps with 
shearless curve at winding number uj (see Refs. Il2ll7ll8h . Operating with 7^[o,i,ii,i,i,...] infinitely many times on the 
standard nontwist map at criticality of the shearless meander studied here limits to a map that is a period- 12 fixed 
point, A, of the renormalization group operator. 

In the vicinity of A, two unstable eigenvalues Si and 82 can be computed to characterize the fixed point. As shown 
in Ref. 0, these are given by 8i = 1 y / TJ17 i , where 



v x = hm -t- 



V2 = lim 



n^oo \ a 



a c[n+12] — O-c 



Or 



(10) 
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FIG. 8: Two levels of magnification, in symmetry line coordinates (x',y'), of the [0, 1, 11, 1,1,1...] meander at breakup. Also 
shown are the nearby up and down orbits of the [12]th (top) and [24] th (bottom) continued fraction convergents. 



where a c \ n ] is the a-value along the bifurcation curve of the [n]th convergent, rather than along the shearless curve, at 
which the sequence of convergent residues exhibits nontrivial limiting behavior (i.e., converges neither to nor oo). 
Using n = 21, we obtain the eigenvalues 

Siks 2.680, <y 2 « 1-584, (11) 

which are, within a small numerical error, the same values that were found for the I/7 2 shearless curve in Ref. 0, and, 
with a slightly larger error, for 1/7 in Ref. HsU£ 



IV. BREAKUP OF THE w = [0, 1, 11, 1, 1, . . .] OUTER SHEARLESS TORI 



A. Search for critical parameter values 



The procedure for finding the critical parameter values for the breakup of outer tori is the same as that described 
in Sec. 1111 A1 The relations y(b), at fixed a, for each of these orbits from Tableware found along symmetry lines. The 
collision parameters 6r n i are again extrema of b{y), only now the outer maxima are used, as illustrated in Fig. [U] for 
the lowest four convergents. The 6[„] of the outer shearless orbits are marked by solid circles (»)^ As in the central 
meander case, the best approximation to 600(a) (here br^i) is used for applying Greene's residue criterion. Finally, 
varying a to find the transition between residue convergence to and 00 results in the parameter values (a co , b co ). 

Even though the 6[ n i(a) seem to follow the same scaling law, Eq. ©, as the ones approximating the central shearless 
meander, the data does not show sufficient evidence for a periodicity of B(n). Using the numerical value for 8\ from 
Eq. |[5J to allow for a direct comparison, and again plotting the offsets of the logarithmic 6[„] differences about this 
average slope, Fig. llOl is obtained»S2, Without an apparent periodicity, a better approximation to 600 than 6[ 34 ] , similar 
to the extrapolation in Eq. ©, cannot be found. 
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FIG. 9: y-values in phase space of the 11/12, 12/13, 23/25, and 35/38 outer orbits on the si symmetry line as a function of 
6, at a = 0.9757564461. 
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FIG. 10: b\ n ] differences in approximating the critical shearless u) = [0, 1, 11, 1,1,.. .] outer orbits at a co = 0.9757564461 on the 
si (above) and S4 (below) symmetry lines. 



B. Residues at breakup 

As before, we search along (a, 6(34] (a)) for the transition between residue convergence to and 00. However, unlike 
before, where the critical parameters could be determined very accurately by observing significant deviations from 
the residue six-cycle for small changes in parameters, here the determination of the exact point of breakup is slightly 
less accurate. Because no clear cycle pattern is apparent, the only criterion to rely on is the residue convergence to 
and 00, which, however, leaves a transition range where the convergence of numerical data is inconclusive. With this 
in mind, our "best guess" for the critical parameters for the outer shearless tori breakup is: 

(aco, b co ) = (0.9757564461, 0.1878717476259388) , (12) 

where the last digit in a co and several of the last digits in b co are merely conjectured to be the values that produce 
the "best transitional behavior" in the residue plots in Fig. ^2 For a better view of the accuracy of these critical 
parameters, all parameters used are listed in Table llVl 

In this case, the only (pairs of) existing periodic orbits [n] of Table [I] are the ones with odd n. Since the outer 
shearless orbits (each considered separately) are not S*-invariant, a residue correspondence scheme as in Table ITT1 docs 
not exist. Therefore residues for both orbits on all four symmetry lines are shown in Fig. 1111 

In contrast to all previous studies of breakups of central shearless tori with noble winding number, for the outer, non- 
S'-invariant tori, no six-cycles occur. Whether the residues converge to higher period cycles remains an open question, 
since within the limited number of 19 data points that could be obtained numerically with sufficient accuracy (up 
to [n] = [37] with only odd [n] existing) no such cycle could be clearly identified, but certainly not ruled out either. 
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a [34] 


6[34] 


0.9757564459 
0.9757564460 
0.9757564461 
0.9757564462 
0.9757564463 


0.1878717467093720 
0.1878717471676554 
0.1878717476259388 
0.1878717480842221 
0.1878717485425054 



TABLE IV: Numerical values of the critical parameters (0(341, 6(34]) in the vicinity of the shearless outer tori breakup, with the 
conjectured critical parameters shown in italics. 



V 


Central torus 


Outer torus 


2 


0.670 ± 0.002 


0.691 ± 0.004 


3 


0.669 ± 0.004 


0.683 ± 0.005 


4 


0.670 ± 0.004 


0.685 ± 0.006 


5 


0.671 ± 0.004 


0.665 ± 0.006 


K x {8) 


0.67 ±0.01 


0.69 ±0.01 


V 


Central torus 


Outer torus 


2 


0.692 ± 0.003 


0.756 ± 0.003 


3 


0.696 ± 0.003 


0.757 ± 0.003 


4 


0.693 ± 0.003 


0.769 ± 0.004 


5 


0.694 ± 0.004 


0.757 ± 0.004 


Ky(0) 


0.69 ±0.01 


0.76 ±0.01 



TABLE V: The regularities of hull functions K x (8) — 8 (top) and K y (8) (bottom), found using the slopes of lines in Fig. 1141 



What can be established, however, is that the outer shearless tori represent the first example where a breakup type 
different from a residue six-cycle is observed in the standard nontwist map. 

C. Phase space at breakup 

The phase space at the critical parameter values for the breakup of the outer shearless tori is shown in Fig. In 
contrast to the case of the central meander no scaling invariance was found. 

V. REGULARITY OF CRITICAL INVARIANT TORI 

The shearless irrational orbits can be described by a function called the "hull function." For an invariant torus of 
rotation number w, this is given by the map K : T — > T x R such that M o K(0) = K(6 + u), and the range of K is 
the invariant torus under consideration. We choose the lift K : R — * R x R of the map K to satisfy 

K(9 + I) = K{9) + (1,0), (13) 

which corresponds to the lift M of the SNM satisfying M(x + l,y) = M(x,y) + (1,0). Equation (|13|) also implies 
that the functions K x (9) — and K y (9) are periodic. These functions for the critical central and outer meanders are 
shown in Fig. 1131 ____ 

We studied these functions using techniques from harmonic analysis developed in Refs. ll9l20l In particular, Fig. 1141 

shows the plot of log \\{d/dt) r 'e- t ^-^I^K(6)\\ L oo versus logi, calculated using 2 25 Fourier coefficients of K{9). 

We see that these functions saturate the bounds given in Ref.lH Ch. 5, Lemma 5: \\(d/dt)^e- t ^- d2 ' d62 K(9)\\ L ^ < 
Ct a ~ v , where a is the Holder exponent of K. This allows us to calculate the regularity of these hull functions from 
the slope of lines in Fig. 1141 and the results are presented in Table Ivl We conclude that the regularity of the central 
shearless torus is 0.68 ± 0.02 while that of the outer torus is 0.72 ± 0.05. 

This agrees very well with the regularity of other shearless noble tori studied in Ref . OS 
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FIG. 11: Residue behavior for orbits on si to S4 at the breakup of u> = [0, 1, 11, 1, 1, . . .] outer shearless tori. Residues for 
assumed breakup parameter values (a co ,b C o) are indicated by the bold line. Residues for five parameter values with a < a co in 
steps of 10 -10 are indicated by thin lines and residues for five parameter values with a > a co in steps of 10 -10 by dotted lines. 
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FIG. 12: Phase space of the [0,1,11,1,1,1...] outer shearless tori at breakup, where a co — 0.9757564461 and b co — 
0.1878717476259388. The symmetry lines si to s 4 are shown. 




VI. CONCLUSION 



In this paper we presented the breakup of two types of shearless invariant tori with noble winding number that 
had not been studied previously: a central meander and an outer torus. The breakup of the central meander showed 
within numerical accuracy the same critical residues, scaling parameters, and eigenvalues of the renormalization group 
operator as the central shearless invariant tori previously studied. From a renormalization group point of view this 
was to be expected: all nontwist maps with a critical shearless torus of noble winding number are expected to be 
equivalent under renormalization to the map with the critical shearless golden mean torus, independent of being a 
meander or not. 

In this light, the result of the outer torus breakup is surprising. Although the winding number is noble, no critical 
residue pattern could be established within the numerically accessible range. This suggests that the number theoretic 
properties of the winding number might not be enough for the classification of different breakup scenarios. In the 
case of nontwist maps the symmetry properties of the shearless torus (here: ^-invariant vs. not S'-invariant) seem 
to affect the breakup. It is possible that after an appropriate coordinate change, that will make the outer torus 
symmetric in those coordinates, the SNM with critical outer torus is equivalent under renormalization to the fixed 
point with critical shearless golden mean torus. Alternatively, this could be the first indication of a new fixed point 
of the renormalization group operator for area-preserving maps. 
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For K y (0) of central torus 
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FIG. 14: The plot of log ||(0/0t)"e 



ri -ty/-d?/de< 



K{Q)\\l<*> versus logi for the hull functions. 
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APPENDIX A: BASIC DEFINITIONS 



For reference, we list a few basic definitions used throughout the main text: 

An orbit of an area-preserving map M is a sequence of points {{xi, yi)}°^ = _ 00 such that M (xi,yi) = (aJi+i, J/i+i)- 
The winding number lu of an orbit is defined as the limit lu = limj_ >00 (xj/z), when it exists. Here the x-coordinate is 
"lifted" from T to R. A periodic orbit of period n is an orbit M n {xi,yi) = (xj +m,yi), Vi, where m is an integer. 
Periodic orbits have rational winding numbers lu = m/n. An invariant torus is a one-dimensional set C, a curve, 
that is invariant under the map, C = M(C). Of particular importance are the invariant tori that are homcomorphic 
to a circle and wind around the x-domain because, in two-dimensional maps, they act as transport barriers. Orbits 
belonging to such a torus gencrically have irrational winding number. 

A map M is called reversible if it can be decomposed as M = I± o I 2 with if = 0. The fixed point sets of h are 
one-dimensional sets, called the symmetry lines of the map. For the SNM the symmetry lines are s\ = {(x, y)\x = 0}, 
s 2 = {(x,y)\x = 1/2}, s 3 = {(x,y)|x = a (l - y 2 ) /2}, and s 4 = {(x,y)|x = a (l - y 2 ) /2 + 1/2}. 

The m/n-bifurcation curve b = $ m / n ,i(a) is the set of (a, b) values for which the m/n up and down periodic orbits 
on the symmetry line Si are at the point of collision. The main property of this curve is that for (a, b) values below 
b = <& m /n,i( a )> t nc r / s periodic orbits with r/s < m/n exist. Thus, m/n is the maximum winding number for 
parameter values along the m/n bifurcation curve. As detailed in Ref. |2, in certain parameter regions multiple orbits 
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of winding number m/n, and therefore multiple bifurcation curves of the same winding number, can exist. 
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